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ABSTRACT 

When  light  propagates  with  the  wave  normal  in  the 
direction  of  an  optic  axis  of  a  biaxial  crystal,  the 
usual  ray  theory  breaks  down.   This  phenomenon  can  be 
analyzed  by  means  of  an  asymptotic  solution  of  Maxwell's 
equations.   The  intensity  is  governed  by  a  partial 
differential  equation  within  the  phase  surfaces,  instead 
of  ordinary  differential  equations  along  rays.   This 
example  shows  that  light  does  not  always  propagate  along 
rays.   A  similar  phenomenon  occurs  in  hydromagnetics. 
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CONICAL  REFRACTION  IK  CRYSTAL  OPTICS 
AND  HYDROMAGNETICS 


1.   Introduction. 

The  phenomenon  of  conical  refraction  is  an  important 
exception  to  the  usual  laws  of  geometrical  optics.   In  the 
case  of  conical  refraction,  a  single  ray,  upon  striking  an 
interface,  splits  into  a  cone  of  rays  (Pigxire  1).   The 
energy  of  the  corresponding  wave  is  usually  thought  of  as 
beiner  transported  along  bundles  of  rays.   When  conical 
refraction  takes  place,  a  bundle  of  rays  splits  into  a 
conical  shell  of  rays  (Figures  2  and  3)«   Nevertheless, 
the  energy  contained  within  the  original  bundle  is  spread 
over  the  entire  region  bounded  by  the  'conical  shell,  not 
merely  the  shell  itself  (as  one  might  believe  from  the 
usual  principles  of  geomietrical  optics). 


Refr^'cted 
rays 


Incoming 
ray 


"^W 


Interface 


f^ 


I 


Bundle 


Figure  1 


Figure  2 
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The  shaded  area  is 
the  locus  of 
endpoints  of 
refracted  rays. 


Figure  3 

This  means  that  the  first  term  of  the  asymptotic 
expansion  of  geometrical  optics  is  different  from  zero 
not  only  in  the  shaded  annulus  shown  In  Figure  3,  but  in 
the  entire  disc  bounded  by  the  outer  circle.   The  mathe- 
matical explanation  for  this  phenomenon  is  that  the 
amplitude  of  the  refracted  wave  Is  not  governed  by 
ordinary  differential  equations  (as  is  usual  in  geometri- 
cal optics),  but  by  a  hyperbolic  partial  differential 
equation  within  the  phase  surfaces  of  the  refracted  wave. 
The  ray  conoid  of  the  partial  differential  equation  is 
identical  with  the  cone  of  refracted  rays  emanating  from 
a  point . 

We  •^hall  demonstrate  the  validity  of  the  preceding 
statements  in  two  cases.   The  first  Is  the  classical  case 
of  propagation  of  light  Incident  on  a  biaxial  crystal  in 
the  direction  of  an  optic  axis.   The  second  case  involves 
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the  propagation  of  a  small  hydromagnetic  discontinuity 
across  a  plane  normal  to  the  magnetic  field,  in  the 
special  case  where  the  Alfven  and  sound  speeds  are  equal. 
In  harmony  with  the  general  rule  that  the  first  term  in 
an  asymptotic  expansion  is  easier  to  compute  than  the 
exact  solution,  each  of  these  problems  can  be  reduced  to 
the  solution  of  the  wave  equation  with  three  independent 
variables. 

An  example,  given  at  the  end  of  the  following  section, 
shows  that  if  the  incoming  wave  vanishes  except  on  a  small 
bundle  of  rays,  then  almost  all  of  the  energy  is  transported 
along  the  refracted  rays. 

It  should  be  pointed  out  that  the  following  analysis 
applies  only  to  propagation  with  the  wave  normal  pointing 
in  a  specific  direction,  i.e.  along  an  optic  axis.   For 
neighboring  directions  of  the  wave  normal,  ordinary 
geometrical  optics  applies.   Thus  there  is  a  change  in  the 
asymptotic  behavior  of  solutions  when  the  wave  normal  is 
near  a  conical  point  of  the  normal  surface.   This  difficulty 
in  the  theory  has  not  been  satisfactorily  resolved,  and 
hinders  comparison  with  experiments. 

2.   Crystal  Optics. 

Consider  the  following  situation:   we  have  a  biaxial 
crystal—'  cut  in  a  plane  normal  to  an  optic  axis.   There 
1/  The  equations  of  crystal  optics  are  discussed  in  [2], 
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is  an  incoming  plane  wave  which  is  normal  to  the  surface 
of  the  crystal,  i.e.  the  incoming  wave  propagates  in  the 
direction  of  the  optic  axis.  Let  the  x-axis  be  an  optic 
axis,  and  the  7,z  plane  be  the  surface  of  the  crystal. 


Incominc  wave 


Figure  l\. 


Let      e,,    e^      and     e^     be   unit   vectors    along   the 
principal   axes    of  the   dieletric    tensor   of   the    crystal, 
and   let      e-,  ,    e-      and      e-^s      be    the    corresponding  eigenvalues, 

We    assume    that      e-i    >    ^p  ^    ^3*      ^^*      ^     ^®    ^^^  magnetic 
permeability.      We   shall  make   use    of  the   following 
abbreviations : 

(i   =   1,2,3)    , 


"i  - 

[I   e^ 

^1  =. 

^3     = 

V  ^3-*l 

1 

Y 


2/57 


V( 62-5^ M 63-62 ) 
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Let      e,  ,    ep      and      e^     be   vmlt   vectors    in  the 
direction   of  the   x,j   and   z   axes,    given  by   the    follovjlng 
equations : 


«1  = 

-  ^^e^    -  0/^ 

^2     = 

t 

"3     = 

-  ^A   ^  ^1^3 

These  equations    fix  the   orientation  of   the   principal 
ajces   of   the    crystal   relative    to   the  x,y  and   z   axes. 

We    asstme    that  the    incoming  wave   has    the    following 
form: 

H     =     H(y,z)    e^*^^*"^^    +   0(i) 
(1)  (x  <    0) 

E     =     I(y,z)    e^'^^^-^^    +   0(1)    . 

Here  the  vectors   H   and   E  represent  the  magnetic 
and  electric  fields;  the  polarization  and  intensity  may 
depend  upon  y   and  z,u>     represents  the  frequency,  and 
the  terms  of  order  l/u)       indicate  that  we  are  interested 
only  in  the  first  term  in  an  asymptotic  expansion  for 
large  oJ    .   The  unit  of  time  is  chosen  so  that  the 
velocity  of  light  in  vacuum  is  unity. 

The  boundary  conditions  on  the  surface  of  the 
crystal  enable  us  to  prescirbe   H   and  E  for  x  =  0. 
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In  order  to   calculate      H      and     E      for     x  >   0,      we    apply 
the   theory   of   asymptotic    solutions   of  partial   differential 
equations    (outlined    in  Section  I4.)    to  Maxv/ell's   equations. 
Accordingly,    we   have 

H     =     H(x,y,z)    exp    iluit ^^)  )    +   0    {h 

(2)  2 

E     =     E(x,y,z)    exp    (i6J(t ^-))    +   0    (i)      , 

for     X  >   0.      H     and     E     must   each  be   a  linear   combination 
of  two  fixed  vectors: 

H(x,y,z)      =      cr^(x,y,z)    H^^^    +   (r^(x,y,z)   H^^^ 
(3) 

E(x,y,z)      =     cr^(x,y,z)   E^^^    +   o-^(x,y,z)   E^^^    . 


Here 


H  -     e^      , 


w 


2      ®2 


f 


H              ~  ®2      * 

(2)  > 

E'     '      =  2  Y   M-   ©1    -   |J- "/el      ©3 

1  2 

The    scalar s      C"         and  (y        must    satisfy  the 


following   system   of  equations: 

-  y  - 


fT~    <JO""  <)0"  _        r\ 


These  equations  are  equivalent  to  the  wave  equation  with 
three  independent  variables.   We  can  see  this  most  easily 
by  introducing  new  variables   4   and   ^  : 

z   =   Y  4  +  ^  . 
Then  we  have 


2 


(4  -  Y:^)  cr^  -  Y<^  =   0 


* 


(6) 


hence 


^  +  ("A  +  Y  4) 


2 


-  ^  ^  "  ^Jt  "  ^  J^^  ^    =   ° 


[-^  -  Y^  (-^  +  ^n  cr^  =  0  , 

Thus  we  see  that  the  first  term  in  the  asymptotic 
expansion  of  the  field  inside  the  crystal  is  detennined 
by  a  solution  of  the  wave  equation  in  three  independent 
variables,  where  the  optic  ajcis  takes  the  place  of  the 
time  axis.   Since  Huygens'  principle  does  not  hold  for 
this  equation,  we  conclude  that  propagation  along  rays 
cannot  take  place.   For  example,  if   H(0,y,z)   and 
E(0,y,z)   vanish  outside  a  small  neighborhood  of  the 
origin,  then  H(x,y,z)   and  E(x,y,z)   will  vanish 
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outside   a  neighborhood   of  a   certain  cone,    and  in  general 
will   be    different  from   zero   inside    the    cone.      The    cone    in 
question   is  the   ray   cone   of    (5),    which  has    the   equation 


X    z      =      z      +7 


The  optic  axis  is  a  generator  of  this  cone,  and  the 
intersection  of  the  plane   x  =  1  with  the  cone  is  a 
circle  of  radius   y/  -/^T  . 


Figure   5 


As  a  (hopefully)  physically  meaningful  example,  we 

1         2 
shall  calculate  (T"         and  cT         in  the  interior  of  the 

cone  of  conical  refraction,  if  there  is  a  small  source 

1         2 
at  the  origin.   The  behavior  of   o~    and   a"    along 

the  refracted  rays  is  more  complicated,  and  depends  upon 

the    shape    of   the   soijrce,    and   the    distribution  of   intensity 
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within  the  source. 

First  we  consider  the  corresponding  problem  for  the 
wave  equation. 

Lemma;   Let   g(a,bj   be  a  continuous  function  such  that 
g(a,b)   =0    if    a^  +  b^  >   1   , 


r(a,b)  da  db  =  1   . 
Let  u(C,>?,^)   satisfy 

u  (o,>^  ,V     =    0 

Then        1.        //  (u^  "•■  u?  +  u|)    d\    d^        is   independent   of 
e      and     £,,      and 


2.      If     Y^   €^    -    (^^   +   ^^)      >      a     >      0    , 


then 


_  1 

271Y 

Proof:   First  we  note  that,  for  £,   =   0    , 

(7)    r/(u|+u^+u^)  dyi    dK     =      /[  g^(a,b)  da  db  . 

Thus  the  energy  of  the  initial  aata  is  independent  of  e. 
But,  for  solutions  of  the  wave  equation,  the  total  energy 
is  independent  of  £,,      This  proves  our  first  assertion. 
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In  order  to   prove   the   second  assertion,    we   represent 
u     in  teiTOs   of   a  fundamental    solution: 


nU,r^,K)     =     -H/( 


g(f,f )    dy   dz 


u(C  yi,V      =       ^       {{ g(a.b)    da   db 

Expanding  the  integrand  with  respect  to  e,  we  obtain  (7). 

Now  we  can  use  our  lemma  to  obtain  approximate  solutions 
of  (5),  or  (6).   First  let  the  data  for  cr'^      and  g-^     have 
the  form 


0-^(0, >|,^)   =   0   , 


where   g(a,b)   is  the  function  to  which  the  lemma  refers, 
and  ^    ~  J'      Then  we  may  set 

2     ^     ^A 


<r 


Y  ^  ^  u,n.,^)  . 


Hence 


(8) 


-'  =   -^^f7^  -0(e2)   , 


r 
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v;here 


r   (x,7,zj      =     ^  xz    -   z2   -  y2 


V^ 


Similarly,    11' 

(r^(0,>^,^J      =     0 


c=r^(o,>^,^j     =    i  g(5,|)     , 


e  °  e'  e 


then 


(9J 


Z-2Y   6'^^^   X 

^     -   2i p;372 —  "  °^'  ^    • 

Combining  18)  and  (9),  we  obtain  a  general  solution  of  (5). 

V7e  may  interpret  these  results  as  follows:   Consider  a 
plane  wave  incident  on  a  crystal  in  the  direction  of  an 
optic  axis,  where  the  amplitude  of  the  wave  vanishes 
except  on  a  bundle  cf  rays  of  size   e.   We  assume  that  e 
is  large  compared  to  the  wavelength  of  the  light,  but 
small  compared  to  the  distance  from  the  interface  to  the 
point  of  observation   (x  »  e  »  l/cJ ) .   If  the  intensity 
of  the  incoming  wave  is  normalized  so  that  the  refracted 
wave  carries  unit  energy,  then  the  amplitude  of  the 
refracted  wave  is  proportional  to   e,   inside  the  cone  of 
conical  refraction. 

On  the  other  hand,  since  the  refracted  wave  carries 
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unit  energy,  we  conclude  that  the  energy  carried  along 
and  near  the  refracted  rays  is   l-0(e  ).   Thus,  although 
the  amplitude  of  the  refracted  wave  is  different  from 
zero  inside  the  cone  of  refracted  rays,  almost  all  of  the 
energy  is  transported  along  rays, 

3,   Hydroir.agnetics . 

In  this  section  we  shall  apply  the  technique  used 
for  crystal  optics  to  the  linearized  equations  of  hydro- 
magnetics,  and  discuss  a  phenomenon  of  conical  refraction 
exactly  analogous  to  conical  refraction  in  crystals. 

We  shall  deal  with  the  equations  for  small  hydro- 
magnetic  disturbances--',  where  the  basic  flow  is  at  rest, 
with  constant  density   p  ,   and  constant  magnetic  field 
vector  H  .   The  equations  for  the  pertiirbations   p,   u 
and  H  are : 

^  ^  p^  div  u  =  0   , 

p„  4^  +  a^  grad  p  -  n  (curl  H)  x  H^  =   0   , 


'o  Jt 


3t 


-  curl  (u  X  H^)   =   0   . 


Here   a  is  the  sound  speed  of  the  underlying  material, 
and  |i,  is  the  magnetic  permeability.   It  will  be 
convenient  to  put  our  equations  in  symmetric  hyperbolic 

—'    These  equations  are  discussed  in  [Ij. 
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form.   Let 


d   =  —  P   , 
Po 


M  =  /-^  H   , 
VPo 

M.  =.M-  H 


o 
Then  our  equations  take  the  symmetric  form 


d.  +  a  dlv  u  =   0   , 

(10)      u.  +  a  grad  d  -  (curl  M)  x  M   =   0 
t  o 


M.  -  curl  (u  X  M  )   =   0 
t  o 


We  look  for  solutions  of  (10)  of  the  form 

(11)      u  I  =  h(t  -  f )  I  u  )+  continuous  function, 

M 

where   h(  )   Is  the  Heavlslde  step  function.   That  Is, 

the  solution  shall  have  a  jump  discontinuity  which 

propagates  with  speed   a   in  the  direction  of  the  x-axls. 

We  make  the  additional  assumption  that   H    (and  hence   M  ) 

o  o 

points  in  the  direction  of  the  x-axis,  and  that  the  Alfven 
speed  is  equal  to  the  sound  speed,  i.e. 


The  propagation  of  such  discontinuities  is  covered  by 
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the  asymptotic  theory  of  partial  differential  equations 
outlined  in  Section  i^..   According  to  the  theory  we  may 
represent  the  discontinuous  part  of  the  solution  in  the 
form 

d' 
a-  (x,y,z,t) 


(12) 


where 


^\ 


/  ji) 


u 


M  / 


/d(l) 


1=1 


u 


\ 


M 


(1) 
(1) 


(2) 


(2) 


u 


M 


(2) 


u 


M 


(3) 
(3) 
(3) 


/  0 


-e 


3 
3 


\ 


u 


M 


(1) 


(1) 


Here,  as  before,   e,  ,  e-   and  e -.      are  unit  vectors 
along  the  x,y  and  z  axes,  respectively.   The  scalars 


cr  ,  (T 


and  cr        must  satisfy  the  following  set  of 


equations : 


-  17  - 


If  we  introduce  new  variables  ^,    x,  y,  t  by  means 
of  the  equations 

t      =     0  +^ 

^        a 

X   =   X 

y  =  y 

z   =   z   , 

then   our   equations   take   the   form 


^21^  +  1  <^o-^ 
2        -,    j_3 


2     dy 


M),  1  ^cr  1  (jo-        _      p^ 

1  <J<r  -'-       1  dcr^    .   jcj-^     =     0     . 

Thus  all  differentiation  takes  place  within  the 

planet,  0  =   const.   It  is  easy  to  see  that  cT    ,    cf 

3 
and  cT -^      satisfy  the  following  equations: 
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J_  (  ^ 


dx     Sx 


T  (^  +  A)) 


o-l      =      0 


cr2     =     0 


s2 


^z 


ax2 


\  (7^  +  -;^))  cr^     =     0 


^1 


H' 


Thus,    as   in  the    case   of   crystal   optics,    the   first 
tenn   in  the   expansion  of  the    solution  is   essentially 
governed  by  the  wave    equation   in  three    independent 
variables.      The   locus    of   a   discontinuity   initially 
localized   at   the    origin  will  be   the    interior   of   a 
circular   cone    in     x,y,z      space. 


Pigijre   6 


-  19   - 


k*      The  Geometrical  Optics  Approximation. 

In  this  section,  we  shall  Indicate  how  the  systems  of 
equations  which  govern  the  first  tenn  of  an  asymptotic 
expansion  are  derived.   In  the  Interests  of  brevity  and 
clarity,  we  consider  a  general  linear  hyperbolic  system 
of  first  order.   This  theory  is  presented  in  greater 
detail  in  [3], 

Too 

Let   u(x  ,x  ,x  ,t)   be  a  vector  with   k  components, 
and  consider  the  linear  system 

(15)         Zu     =   u,  -K  ^  A^  ^  +  Bu  =   0   . 

v=l     ^X^ 

Here   A    and  B   are   k  x  k  matrices.   In  the  case  of 
crystal  optics,   u  represents  the  pair  of  vectors   (E,H) 
and  oCu  =  0   corresponds  to  Maxwell's  equations.   In  the 
hydromagnetic  case,   u  represents  the  triple   (d,u,M) 
and  r>C  V.  =   0      corresponds  to  the  system  (10).   We  look 
for  formal  solutions  of  (IS)    of  the  form 

(16)         u  =  ^«l^'-.t' 


III 1 —  aJ(x,t) 

J=0    (iU)J 


For   large   values    of  the   frequency      cj    ,      the    first   term 
in  the    expansion   is   a   good   approximation  to   the    solution 
u;      we   vfrlte 


u      ^J 


^lujflU.t)  ^o,^^^,      _ 
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Upon  substitution  of  (16)  into  (15),  we  see  that 
the  following  set  of  equations  must  be  satisfied: 

(17)  A  a°   =   0   , 

(18)  A   a^""^   +    Xuh      =      0  (j   =   0,1,...)    . 
Here      A      is    the    characteristic  matrix: 

Equation  (1?)  implies  that  the  determinant  of  A  must 
vanish;  hence  0(x,t)  must  satisfy  the  characteristic 
equation. 

We   assume    that,    for  fixed     ^,      there    are   exactly 
r      linearly   independent    solutions    of   equation    (17): 

A  R^^^      =0  (i   =  l,...,r}    . 

Then  we  must  have 


(19)  a°(x,t)   =  tZ   <^^(x,t)  R^^^   , 

i=l 

where  cr        are  scalar  factors.   There  are  also  r 
linearly  independent  vectors  L     such  that 

L^^^  A  =   0  (1  =  l,...,r)  . 

Now  consider  equation  (18),  setting   J  =  0.   A  necessary 
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and  sufficient  condition  for  this  equation  to  have  a 
solution  is  that 


l'^'  =^(a°)  =  0        (j=i r)    , 


or 


(20)      l(3'^(  tl'y-'u.t)   r'I',  =0   (J  =l,...,r)  . 

The  system  (20)  corresponds  to  the  system  (5)  in  the  case 
of  crystal  optics,  and  the  system  (13)  in  the  case  of 
hydromagnetics.   Thus  the  leading  term  in  the  expansion 
of  u   is  determined  by  (20),  together  with  appropriate 
initial  and  boundary  conditions. 

It  has  been  customary,  both  in  the  mathematical  and 
physical  literature,  to  treat  only  those  cases  in  which 
(20)  reduces  to  a  system  of  ordinary  differential 
equations,  or  a  single  equation.   In  such  cases  the 
unknowns  c         are  differentiated  along  rays  and  the 
energy  is  indeed  transported  along  bundles  of  rays. 

When  conical  refraction  takes  place,  the  system  (20) 
is  a  nontrivial  hyperbolic  system,  and  the  cone  of  rays 
represented  in  Figure  1  corresponds  to  the  ray  conoid  of 
system  (20).   This  fact  can  be  verified  directly  in  the 
cases  of  crystal  optics  and  hydromagnetics.   The  general 
theorem  follows  from  a  lemma  of  P.  D.  Lax  and  R.  Courant, 
quoted  in  [3 J,  pp.  33-3k* 
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